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1.  OBJECTIVE  uui s^ibinrEi 

'The  objective  of  this  MTP  is  to  determine  the  effect  of  high  temp^-a- 


o 

*^ature  and  temperature  changes,  resulting  from  aerodynamic  heating,  on  the  missile 
v^^skin,  structure,  and  interior  components  using  laboratory  methods.  Such  inform- 
\  nation  is  used  to  determine  the  adequacy  of  equipment  to  meet  specified  require- 
‘■‘•srments  such  as  the  ability  to  withstand  the  aerodynamic  heating  associated  with 
00  flying  a  specified  trajectory. , 

■J  2.  BACKGROUND 


Q 


Aerodynamic  heating  is  the  conversion  of  kinetic  energy  into  heat 
energy  as  the  result  of  the  relative  motion  between  a  body  and  a  fluid,  and 
the  subsequent  transfer  of  this  heat  energy  through  the  skin  into  the  structure 
and  the  interior  of  the  vehicle.  Some  heat  is  produced  by  fluid  compression  at 
and  near  stagnation  points,  such  as  the  missile  nose  or  the  airfoil  leading 
edges.  Additional  heat  results  from  friction  along  the  missile  skin  inside  the 
boundary  layer.  In  this  area,  the  relative  motion  between  the  fluid  particles 
and  the  vehicle  is  reduced  by  shear  force. 


Since  the  advent  of  high  speed  vehicles,  the  effects  of  aerodynamic 
heating  have  been  of  extreme  importance  to  the  design  engineer.  Frictional 
heating  becomes  significant  at  velocities  of  approximately  2000  feet  per  second 
(fps),  which  are  well  within  the  capabilities  of  some  of  the  current  missiles. 

It  increases  the  problem  of  mechanical  load  distribution  of  the  vehicle,  weak¬ 
ening  structural  members  so  that  more  weight  is  required  to  withstand  a  given 
stress.  In  addition,  thermal  stresses  further  modify  the  mechanical  stress 
pattern  because  of  the  differential  in  expansion  rates  for  dissimilar  materials. 

The  test  engineer  must  simulate,  in  the  laboratory,  many  test  condi¬ 
tions  prescribed  either  by  the  designer  of  the  equipment  or  the  military  sped- 
fications . 

3.  REQUIRED  EQUIPMENT 

a.  Applicable  Testing  Facility  (convective  or  nonconvective ) 

b.  Heat  Control  Computer 

c.  Applicable  Instrumentation  (high  temperature  strain  gauges, 
thermocouples,  pyrometers,  etc.  (as  required)) 

d.  Applicable  Data  Recording  and  Processing  Equipment 
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SCOPE 


5.1  SUMMARY 


This  MTP  describes  methods  of  subjecting  a  test  specimen  to  heating 
effects  that  simulate  those  aerodynamic  heating  effects  that  the  test  specimen 
would  encounter  if  flown  in  a  given  trajectory. 


Included  in  the  MTP  are  methods  for  mathematically  determining  the 
probable  heating  effects  on  a  test  specimen  flying  a  given  trajectory  using 
standard  air  tables,  the  known  trajectory,  the  shape  of  the  test  specimen  and 
known  heat  transfer  constants. 


Also  included  are  methods  for  simulating  the  above  heating  effects  as 
a  function  of  time  using  existing  facilities  so  that  the  effects  on  the  test 
specimen  from  such  exposure  can  be  ascertained. 

5.2  LIMITATIONS  < 

]>ie  to  the  variety  of  missile  configurations  and  trajectory-velocity 
conditions,  no  attempt  is  made  in  this  MTP  to  describe  a  testing  procedure  for 
any  particular  missile.  The  considerations  outlined  may  be  adapted  as  necessary 
to  apply  to  a  given  missile  and/or  flight  conditions. 


The  mathematical  considerations  presented  in  the  calculation  of  pre¬ 
dicted  heating  effects,  limit  this  test  to  simulated  temperature  and  altitudes 
where  conventional  gas  dynamics  solutions  are  valid  or  which  are  in  the  regime 
of  continuum  flow  (assumed  up  to  an  altitude  of  approximately  30  to  4o  miles) 
and  where  dissociation  and  ionization  are  negligible.  The  test  is  also  limit¬ 
ed  to  applications  where  it  can  be  assumed  that  there  is  locally  (0)  constant 


6.  PROCEDURES 


6.1  PREPARATION  FOR  TEST 
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Personnel  involved  in  testing  must  be  familiar  with  the  missile 
configurations  they  are  to  test,  and  fully  qualified  in  the  use  of  the  assoc¬ 
iated  test  equipment  and  facilities.  Prior  to  the  test,  pertinent  technical 
manuals,  manufacturer's  manuals,  or  other  material  available  must  be  reviewed 
to  make  a  proper  selection  of  test  equipment  and  accessories. 

The  references  listed  in  this  MTP  contain  more  detailed  information 
if  it  is  required. 

6.2  TEST  PROCEDURES 

a.  Determination  of  mission  profile. 

1)  Calculate  the  velocity  and  altitude  of  the  test  specimen  as  a 
function  of  time  for  the  specified  trajectory. 

2)  Using  standard  atmosphere  tables  (see  reference  Y)  determine 
the  following  free  stream  parameters  for  the  altitudes  deter¬ 
mined  above:  Tea,  p«,  Pooi  and  from  the  velocities: 

(See  Appendix  A,  page  A-l). 

3)  Using  the  shape  of  the  test  specimen,  its  attitude  in  flight, 
and  the  known  free  stream  parameters,  determine  the  potential 
flow  parameters,  T0,  V0,  p~,  p0,  and  Mq  using  the  references 
given  in  Appendix  A,  page  A-2. 

4)  Calculate  the  local  heat  transfer  function  (h)  and  recovery 
temperature  (Tr)  (as  defined  in  Appendix  A,  page  A-3 )  as 
functions  of  time  using  the  potential  flow  parameters,  heat 
transfer  constants,  and  the  test  specimen  characteristic 
length.  The  material  and  references  listed  in  Appendix  A 
shall  be  used  to  accomplish  this. 

b.  Determine  the  type  of  testing  facility  to  be  used  in  conducting 
the  test  in  accordance  with  Appendix  B. 

NOTE:  Ensure  that  the  capacity  of  the  testing  facility  is  sufficient 
to  supply  the  maximum  heat  flow  required. 

c.  Determine  the  types  of  instrumentation  and  data  processing  equip¬ 
ment  to  be  used  in  accordance  with  Appendix  C. 

d.  Load  the  heat  control  computer  memory  with  the  values  of  h  and  Tr 
in  accordance  with  Appendix  D  so  that  the  net  heat  flow  into  the  test  specimen 
as  a  function  of  time  will  be  the  same  as  the  test  specimen  would  experience  in 
flying  the  selected  trajectory. 

e.  Verify,  using  a  data  printout  or  other  applicable  means,  that  the 
values  of  data  stored  in  the  computer  memory  by  the  above  step  corresponds  to 
the  pre -calculated  values. 

f.  Place  the  test  specimen  into  the  test  facility. 

g.  Instrument  the  test  specimen  to  monitor  and  record  wall  tempera¬ 
ture  (Tw),  and  other  desired  measurements,  e.g.,  surface  strain  (See  Appendix  C) 
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h.  Operate  the  test  facility  using  the  heat  control  computer  to 
regulate  the  flovT  of  heat  to  the  test  specimen. 

NOTE:  The  onboard  equipment  shall  be  operated  during  the  test  and 
any  variation  from  design  capabilities  shall  be  noted. 

i.  Using  the  test  facility  instrumentation  instead  of  the  test  speci¬ 
men  instrumentation  monitor  the  net  heat  flow  into  the  test  specimen. 

j.  Monitor  all  test  specimen  instrumentation  data. 

k.  Shut  down  the  test  facility. 

l.  Remove  the  test  specimen  and  inspect  for  any  effect  of  the 
exposure  to  the  test  environment.  Record  any  observations,  e.g.,  degradation 
of  material  properties,  creep,  fatigue,  thermal  stress. 

m.  Retain  all  data  recordings. 

n.  The  onboard  equipment  shall  be  operated  and  any  variations  from 
design  capabilities  shall  be  noted. 

6.3  TEST  DATA. 

Retain  all  data  recordings. 

Retain  all  physical  observations  made  on  the  test  specimen  after 
completion  of  the  test. 

Retain  all  observations  of  the  operation  of  onboard  equipment  during 
and  after  the  test. 

6.4  DATA  REDUCTION  AND  PRESENTATION 

Compare  the  resultant  heating  effect  with  expected  values  to  ensure 
that  test  conditions  are  representative  of  calculated  flight  conditions. 

From  observations  made  during  the  physical  inspection  of  the  test 
specimen  and  from  observation  made  of  the  operation  of  onboard  equipment  during 
and  after  the  test  determine  if  the  test  specimen  meets  manufacturer's  specifi¬ 
cations  and  if  it  is  suitable  for  its  intended  usage. 
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GLOSSARY 

1.  Body  of  Revolution:  A  shape  generated  by  a  straight  or  curved  line  turn¬ 
about  an  axis. 

* 

2.  Constant  Pro-pert y  Values:  An  equation  of  property  values  (heat  transfer 
rate,  absolute  temperature,  pressure,  etc.)  considered  to  be  constant  for 
purposes  of  calculation. 

3.  Continuum  Flow:  Particle  movement  in  a  homogeneous  fluid. 

4.  Conventional  Gas  Dynamics:  A  study  of  fluid  in  motion. 

5.  Dynamical  Variables:  Those  variables  in  terms  of  which,  classical  mechan¬ 
ics  are  built,  and  which  can  be  given  an  operational  definition. 

6.  Enthalpy  (H):  H  =  U  +  p V,  where  U  is  the  internal  energy  of  the  system,  p 
is  the  pressure,  and  V  is  the  volume.  At  constant  pressure,  the  change  in 
enthalpy  is  equal  to  the  heat  flow  into  the  system  from  the  surroundings. 

7.  Local  Heat  Transfer  Coefficient  (h):  h  -  British  Thermal  Units  divided 
by  the  product  of  time  in  seconds,  area  in  square  feet,  and  temperature 
in  degrees  Rankine. 

8.  Newtonian  Flow  Theory:  An  expression  of  shear  forces  in  viscous  layers, 
shock  detachment  parameters,  and  boundary  layer  characteristics. 

9.  Nusselt  Number:  The  nondimensional  parameter,  defined  as  Nu  =  £ 

AT  k 

where  Q  is  the  rate  of  heat  loss  from  a  body,  AT  is  the  difference  of 
temperature  between  the  body  and  its  surroundings,  d  is  the  scale  size  of 
the  body,  k  is  the  thermal  conductivity  of  the  surrounding  fluid, 

10.  Prandtl  Number:  The  kinematic  viscosity  of  a  fluid  divided  by  its  thermal 

conductivity,  Pr  =  SBi. 

k 

11.  Recovery  Factor:  A  quantity  used  to  specify  recovery  temperature,  and  how 
closely  recovery  temperature  approaches  free  stream  stagnation  temperature. 

12.  Recovery  Temperature:  In  the  case  of  an  insulated  plate,  the  equilibrium 
temperature  at  which  heat  transferred  out  from  the  inner  viscous  layer  is 
in  balance  with  the  viscous  work  done  on  the  layers. 

13.  Reynolds  Number:  The  ratio  of  the  inertia  forces  to  the  viscous  forces  in 
a  fluid,  Re=  EXi. 

P- 

14.  Runge-Kutta  Method:  A  numerical  method  for  solving  a  differential  equa¬ 
tion. 

15.  Slip  Flow:  One  of  four  distinct  flow  regimes.  In  order  of  increasing 
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free  mean  paths,  they  are  continuum  flow,  slip  flow,  intermediate  flow, 
and  free  molecule  flow.  Information  on  skin  friction  and  heat  transfer 
in  slip  flow  is  limited,  and  hardly  any  theoretical  results  are  available 
for  the  intermediate  regime  between  slip  flow  and  free  molecule  flow. 

16.  Specific  Heat;  The  quantity  of  heat  required  to  raise  the  temperature  of 
a  unit  mass  one  degree. 

17.  Stagnation  Point:  A  point  at  which  the  velocity  of  a  fluid  is  zero  with 
respect  to  a  body  in  motion. 

18.  State  Variables:  A  limited  set  of  dynamical  variables  of  a  system, 
sufficient  to  specify  the  state  of  the  system  for  given  considerations. 

19.  Static  Temperature  (wall  and  free  stream):  An  unchanging  temperature 
which  is  unaffected  by  heat  from  viscous  shear. 

20.  Taylor-Maccoll  Equation:  An  equation  to  determine  the  air  pressure  on  the 
surfaces  of  cones  under  conditions  of  supersonic  flow  at  zero  angle  of 
attack. 

21.  Transition  Point:  The  location  on  a  moving  object  at  which  the  surround¬ 
ing  fluid  ceases  to  follow  the  surface,  developing  eddies  and  random 
motion. 
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APPENDIX  A 

DETERMINATION  OP  THERMAL  ENVIRONMENT  PROM  A  MISSION  PROFILE 

In  the  attempt  to  solve  flow  and  heat  transfer  problem  for  bodies 
moving  at  supersonic  speeds,  the  concept  of  boundary  layers  has  proved  to  be 
of  great  practical  value .  It  permits  treatment  of  extremely  complicated 
actual  physical  flow  fields  in  two  separate  parts: 

a.  Potential  flow  outside  the  boundary  layer  where  shear  forces  in 
the  fluid  can  be  neglected,  compared  to  inertial  forces. 

b.  Boundary  layer  flow,  where  shear  as  well  as  inertia  stress  must 
be  considered. 


Differential  equations  describing  real  physical  flow  have  not  yet 
been  solved.  Applied  to  the  potential  flow  fields,  however,  these  equations 
can  be  greatly  simplified,  since  shear  forces  are  ignored.  The  governing 
equations  of  the  boundary  layer  are  also  subject  to  simplifications,  because 
of  a  simple  flow  geometry. 


n 


All  relationships  described  in  this  section  of  the  MTP  pertain  to 
conventional  gas  dynamics  or  the  regime  of  contiuum  flew,  assumed  to  be  valid 
up  to  an  altitude  of  approximately  30  to  40  miles.  The  correlations  hold  true 
for  bodies  flying  under  zero  or  small  angles  of  attack,  with  locally 
constant  temperature  along  the  surface,  and  with  zero  pressure  gradients  per¬ 
pendicular  to  the  missile  skin,  i.e.,  =  0.  Neither  problems  in  slip  flow 

(free  molecule  flow)  nor  the  temperature  regions  where  dissociation  and 
ionization  occur  will  be  discussed. 


The  prerequisite  for  engineering  calculations  of  temperatures  or  heat 
fluxes  resulting  from. aerodynamic  heating  is  the  knowledge  of  the  trajectory  of 
the  vehicle  to  be  investigated  for  any  particular  mission.  Prom  a  given  flight 
path,  the  altitude  and  velocity  of  the  missile  can  be  established  as  a  function 
of  time  for  any  desired  time  interval.  The  required  free  stream  (°°)  parameters 
are  extracted  from  tables  of  the  Standard  Atmosphere.  (See  Reference  Y).  The 
parameters  are: 


Va»;  Vehicle  Velocity  in  feet  per  second. 

T^j  Atmospheric  Temperature  in  degrees  Rankine 
p^j  Atmospheric  Density  in  pounds  per  cubic  foot. 
Atmospheric  Pressure  in  pounds  per  square  foot. 

Vai 

The  Mach  number  can  be  obtained  as:  M»  =  ___ 
where:  am=  ^MgRT 

QfoB=  4-9.02  (velocity  of  sound  in  fps) 

The  constant  49.02  renders  this  equation  valid  for  air  only. 

If  p  is  not  available,  it  can  be  determined  by: 

P®  *  (53.3T  x  T* 


(A-l) 


(A-2) 


(A-3) 
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(a)  Potential  Flow  Solutions  —  The  next  step  is  to  find  the  para¬ 
meter  values  just  outside  the  boundary  layer  and  behind  the  shock  front 
occurring  on  a  body  traveling  at  supersonic  speeds.  These  values  can  be 
obtained  by  solving  the  potential  flow  equations  for  the  body  shape  to  be 
investigated.  Information  on  potential  flow  solutions  may  be  found  in 
References  A,  B,  C,  and  D.  This  literature  gives  solutions  of  potential  flow 
around  cones,  ogives  and  arbitrarily  shaped  bodies  of  revolution  from  low  to 
high  Mach  numbers.  Reference  A  discusses  the  potential  flow  in  form  of  the 
Taylor-Maccoll  equation  and  its  solution  in  series  by  iteration.  It  also 
explains  in  detail  the  methods  of  characteristics  for  an  axially  symmetric 
supersonic  flow. 

Reference  B  contains  the  calculation  of  second-order  supersonic  flow 
past  bodies  of  revolution  at  zero  angle  of  attack.  The  solutions  are  valid  for 
Mach  numbers  up  to  approximately  five  and  for  arbitrarily  shaped  bodies  of 
revolution. 

The  solution  of  supersonic  flow  fields  by  means  of  the  conical  shock 
expansion  method  above  Mach  five  is  discussed  in  Reference  C.  Reference  B 
presents  a  collection  of  flow  theories  including  flow  over  blunt  nosed  bodies 
at  high  velocities.  The  solutions  for  potential  flow  enable  determination  of 
the  following  flow  parameters  for  the  station  on  the  missile  for  which  a 
temperature-time  history  must  be  established: 

Tq  ;  Jkj  Po  ;  Mq  ;  j3o_ 

T«b  Vg,  Poo  Moo  B® 

With  flow  values  calculated  Just  outside  the  boundary  layer  (o),  the  Reynolds 
number  can  be  calculated  as  follows: 
p  -V  «1 

Re  .  -2 — 2 —  (A-4) 

1*0 

where: 

1  =  the  characteristic  length  in  feet 

p,  =  the  dynamic  viscosity  and  a  function  of  temperature 
expressed  by  a  power  law  given  in  Reference  E: 

0.1458  x  (4  iy*6  , 

“  ■  0.6TS9  x  10-  n0-t'  .  fr  T)  (A-5) 

Next,  determine  whether  the  boundary  layer  is  laminar  or  turbulent  at 
the  point  of  investigation.  At  the  present,  predictions  of  the  transition  point 
location  must  be  based  mainly  upon  empirical  information.  Detailed  studies  are 
presented  in  References  F,  G,  H,  and  I.  As  a  conservative  estimate  and  for 
engineering  calculations ,  however,  assume  that  turbulent  flow  exists  above  a 
local  Reynolds  number  of  1  x  10°  to  2  x  10  . 

(b)  Basic  Relationship  for  Energy  Balances  and  Heat  Transfer  — 
Reference  J  outlines  the  heat  energy  balance  and  evolving  skin  temperature 
equation  as  follows. 


A-2 
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The  local  heat  transfer  rate  into  or  out  of  a  surface  element  of  the 
vehicle  is  equal  to: 


‘local 


=  %  +  *E  +  %  “ 


where : 


q  =  heat  transfer  rate  [2Si] 

sec 

-  heat  transferred  to  the  skin  through  the  effect  of 
aerodynamic  heating 

qE  =  heat  generated  by  equipment  inside  the  vehicle 
qN  =  solar  and  nocturnal  radiation 

q^  «  heat  transferred  from  the  boundary  layer  and  skin  by  radiation 

For  the  domain  of  operation  of  the  vehicle  to  be  discussed,  qN  and  q_  generally 
are  negligible  compared  to  the  other  two  terms,  and  will  be  omitted  In  future 
calculations.  Heat  dissipation  by  radiation  from  the  skin  is  effective  only 
at  relatively  high  temperatures,  as  shown  later  in  the  radiation  term.  The 
equation  for  4iocal  can  be  written  as: 


local 


=  qv  ‘ 


(A-7) 

(A-8) 


where: 


q.  .  *  h  (T  -  T  )  -  soT  4  (A-S) 

*local  '  r  w'  w  v  ' 

h  =  local  heat  transfer  coefficient  [ - — ] 

sec  ft  R 

e  =  surface  emissivity  -  dimensionless 
a  =  Stephan -Boltzmann  constant  C'"ec'  '^a— ] 

Tw  *  wall  temperature  [°R] 

Tr  =  temperature  of  the  gas  in  the  boundary  layer  on  the  skin  surface 
[°R],  when  the  convective  heat  transfer  at  this  point  is  zero. 


Conversely, 


‘local 


=  c-s-6  (£5i) 
dt 


(A-9) 


where: 

BTU 

c  =  specific  heat  of  the  skin  material,  r  ] 

s  =  skin  thickness  in  feet 

6  =  specific  weight  of  the  skin  material  in  pounds  per  cubic  foot 
Ty  =  Wall  temperature,  [°R] 

The  assumed  conditions  are  that  the  inside  of  the  skin  is  insulated  and  that 
all  heat  energy  transferred  to  the  skin  from  the  boundary  layer  is  used  to 
raise  the  enthalpy  of  the  skin,  and  that  the  skin  is  thin  and  Tw  exists 
throughout  the  skin.  These  conditions  do  not  exist  in  reality,  since  part  of 
the  heat  is  transferred  into  the  interior  of  the  vehicle  and  the  temperatures 
obtained  under  those  restrictions  are  conservative.  The  general  expression  for 
the  local  temperature  change  of  the  skin  with  respect  to  time  finally  reads: 


where: 


121*1 1  (T 
dt  G  s  r 


(A- 10) 
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BTU 

G  =  SC6=  heat  absorption  capacity  of  the  skin 

e  =  surface  emissivity  [dimensionless] 

a  =  Stephan -Boltzmann  radiation  constant  and  is  equal  to 
17.3  x  10"10  r  BTU 
°  ~  3^00  '  '-sec  ft*  °R  ■* 

=  wall  temperature  [°R] 

Tr  =  temperature  of  the  gas  in  the  boundary  layer  on  the  skin  surface, 
[°R]  when  the  convective  heat  transfer  at  this  point  is  zero 

Equation  (A-IO)  is  the  basic  relationship  for  the  change  of  skin  temperatures 
with  respect  to  time.  To  this  point,  the  preparatory  work  holds  true  for 
laminar  as  well  as  for  turbulent  boundary  layers.  From  now  on,  however, 
laminar,  turbulent,  and  stagnation  point  heat  transfers  must  be  discussed 
separately. 

(c)  Heat  Transfer  in  Laminar  Boundary  Layers  —  The  main  problem  of 
solving  equation  (A-10)  is  to  determine  the  local  heat  transfer  coefficient  h. 
In  the  case  of  the  laminar  boundary  layer,  numerous  solutions  of  the  differ¬ 
ential  equations  describe  flow  and  energy  conditions  in  the  continuum  flow 
regime.  Three  of  these  are  presented  here  from  Reference  K.  They  are  restrict 
ed  to  two-dimensional  flow  over  surfaces  along  which  the  pressure  is  constant 
in  flow  direction.  References  L  and  M  outline  procedures  to  convert  informa¬ 
tion  obtained  for  two-dimensional  flow  to  conditions  on  cones.  According  to 
these  references  the  local  heat  transfer  coefficient  on  a  flat  plate  must  be 
multiplied  by  to  obtain  the  local  heat  transfer  on  a  cone  at  zero  angle  of 
attack  for  laminar  flow.  For  turbulent  boundary  layers  the  local  heat  trans¬ 
fer  coefficient  on  a  cone  is  found  by  dividing  the  cone  Reynolds  number  by  2 
and  determining  the  local  heat  transfer  coefficient  for  the  flat  plate.  This 
then,  is  the  value  for  the  cone  heat  transfer  coefficient. 

(d)  Solution  for  Constant  Property  Values  --  For  vehicles  traveling 
in  the  transonic  and  low  supersonic  range,  this  method  is  the  simplest  since 
the  property  values  (viscosity,  heat  conductivity,  density,  and  specific  heat) 
are  assumed  to  be  constant.  As  long  as  the  ratio  of  wall  to  free  stream  static 
temperature  lies  between  0.5  and  1.5  the  results  obtained  with  this  method  are 
sufficiently  accurate. 


where: 


This  temperature  can  be  obtained  through  the  following  relation: 


S  S  i 


(a- 11) 


r  =  recovery  ^factor  -  see  Glossary  item  11 
g  =  32.2  Cg'ecg""]  (gravitational  acceleration) 

T0  =  local  temperature  just  outside  the  boundary  layer 

J  =  778.26  [  (mechanical  equivalent  of  heat) 

T^0=  free  stream  stagnation  temperature 


A-^ 
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Solving  for  Tp  yields: 

Tr  =  Tc  +  r  . 


50120  x  c. 


(A -12) 


The  recovery  factor  r  can  be  approximated  as  the  square  root  of  the  Prandtl 
number,,  so  that: 

T  =  T  +  -\Tp  V  a 

r  0  ^  ***  _ 12 _  (a-13  ) 

50120  x  cp  1  i} 

The  local  heat  transfer  coefficient  h  can  be  expressed  in  any  of  the  following 
ways  by  nondimensional  numbers: 

Nu  =  lA~  (k-lk) 


St  = 


p  cp  Vo 


(A-15) 


St  = 


where: 


2  x  P_ 


St  - 


R  P 
e  r 


h  =  convective  heat  transfer  coefficient,  [ - —21 - ] 

*  uhr.  ft*  °F 

k  =  thermal  conductivity  coefficient,  [ — SSI - ] 

hr.  ft 


(A-16) 

(A-17) 


C-n  =  specific  heat  [-—Siii— — 

*  lb .  P  * 

1  =  characteristic  length  in  feet 
p  =  density  [^-  ] 

V0  =  velocity  [-£M 
sec 

Re  =  Reynolds  number 

p.  =  Coefficient  of  viscosity 

Pr,  Cf,  and  ^  are  described  in  equations  A-l8,  A-20,  and  A-21. 

The  Prandtl  number  is  a  property  value  which  relates  viscosity,  specific  heat, 
and  thermal  conductivity  as  follows: 

^  (A-18 

As  mentioned  earlier,  the  recovery  factor  r  can  be  approximated  as: 

r  =  (A- 19 

for  laminar  flow  in  the  range  of  Prandtl  numbers  between  0.5  and  5. 

The  local  skin  friction  factor  is 
0.664 

f  ~  (A-20 

and  the  local  Nusselt  number  is _ 

Nu  =  0.332  x  V*R^  x  €  P^  (A-21' 


pr  = 


Cf 


(A-18) 


(A-19) 


(A-20) 

(A-21) 
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Substitutions  with  the  parameters  discussed  enable  derivation  of  the  expression 
for  the  local  heat  transfer  coefficient: 


h  =  0.332  c. 


x  pr 


2/3 


and  finally: 
dTw 

—  =  0.332 
dt 


P  vc 


"273 


Cr  r 


T  + 
w 


nr 

V 


6  O  Tt 


50120  X  c 


(A-22) 


(A-23) 


P 


This  equation  can  be  integrated  by  the  Runge-Kutta  method,  since  all  parameters 
(once  chosen  for  the  corresponding  temperature)  are  constants  except  the 
velocity.  It  simplifies  numerical  integration  if  the  velocity  can  be  expressed 
analytically  as  a  function  of  time. 


The  advantage  of  this  method  is  that  it  depends  upon  a  minimum  of 
parameters  and  applies  in  general  to  any  fluid. 


(e)  Solutions  for  Proportionality  Between  M-,  k,  and  Tw  —  At  higher 

supersonic  velocities,  the  temperature  gradients  in  the  boundary  layer  perpen¬ 
dicular  to  the  surface  are  so  large,  the  temperature  dependence  of  the  property 
values  must  be  considered.  Heat  conductivity  and  viscosity  vary  more  rapidly 
than  Prandtl  number  and  specific  heat.  Assuming  that  Pr  and  cp  are  constant, 
the  following  expression  shows  that  the  viscosity  varies  proportionally  to  the 


conductivity. 

JL  =  Jk. 

k  cp 


where: 

k  =  thermal  conductivity. 

Reference  N  shows  that  the  relations  for  constant  property  values  can  be  used, 
if  the  temperature  dependent  property  values  are  introduced  at  a  reference 
temperature,  which  is  expressed  as  follows: 

-ft  =  1  +  0.032  x  +  0.58  (SL.  -  1)  (A-24) 

Ao  ^0 

This  enables  performance  of  the  same  operations  Indicated  in  the  preceding 
paragraph. 


f .  Solutions  for  Air  with  Actual  Temperature  Variation  of  all 
Properties  —  The  previous  considerations  were  valid  for  moderately  high 
temperatures  iri  the  boundary  layer,  and  it  was  assumed  that  the  Prandtl  number 
and  the  specific  heat  for  air  were  constant.  This  statement  does  not  hold  true 
for  extremely  high  air  velocities  where  the  Prandtl  number  and  specific  heat 
must  be  considered  to  be  variable  also. 

Better  results  are  obtained  in  the  calculation  of  the  heat  transfer 
coefficient  when  the  property  values  and  state  variables  are  based  on  enthal¬ 
pies  instead  of  temperatures.  The  use  of  enthalpies  is  therefore  recommended 
whenever  the  temperature  variation  in  the  boundary  layer  is  large. 


A-6 
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For  variable  properties,  the  relationship  (A-20)  is  still  valid;  i.e., 
the  friction  coefficient  varies  proportionally  vith  the  inverse  of  the  square- 
root  of  the  Reynolds  number.  However,  it  is  necessary  now  to  .introduce  the 
property  and  state  values  at  a  reference  enthalpy  (i),  as  defined  in  the 
following  form: 

1*  -  i0  +  °-5°  Uw  -  *o>  +  °'22  (*r  -  io)  (A-25) 

Results  obtained  by  using  this  enthalpy  formula  closely  agree  with  the  exact 
boundary  layer  solutions.  In  the  Reynolds  number  g,  and  p  must  be  introduced 
at  1*,  and  the  velocity  at  V0.  The  recovery  factor  is  expressed  by  equation 
(A- 19)  when  the  Prandtl  number  is  based  on  the  reference  enthalpy  (A-25). 


Values  for  the  Prandtl  number,  viscosity,  and  specific  heat  can  be 
found  in  gas  tables  (References  E,  0,  and  P).  The  difference  in  Prandtl  numbers 
found  in  reference  literature  is  considerable.  Use  the  tables  in  Reference  E. 
Using  enthalpies  instead  of  temperatures,  equation  (A-23)  now  can  be  modified 
to: 


0.332  x  p  xV  (i  -  i  +^/T^,V  a  )  -  e  a  T  4 

_  rT  o  o  w  r*  o  _ w_ 

•t  'o  *  VTC* «  (Pr»)2/3  50120  G 


(A-26) 


The  density  pr  is  introduced  at  recovery  enthalpy,  while  all  para¬ 
meters  with  the  superscript*  are  based  oh  the  reference  enthalpy  (A-25).  If 
the  specific  heat  c  of  the  skin  material  varies  considerably  with  temperature, 
the  change  should  be  taken  into  account  and  all  values  for  c  should  be  chosen 
at  wall  temperature  Tw. 

g.  Heat  Transfer  in  Turbulent  Boundary  Layers  —  In  contrast  to  the 
treatment  of  laminar  boundary  layers,  the  calculation  of  friction  factors  and 
heat  transfer  coefficient  in  the  turbulent  boundary  layer  is  extremely  diffi¬ 
cult,  because  knowledge  of  the  turbulent  exchange  mechanism  under  widely 
varying  temperatures  is  limited  and  experiments  must  be  relied  upon. 

A  comparison  between  predictions  of  different  theories  has  shown 
large  differences  and  therefore  a  calculation  procedure  is  presented  which  is 
recommended  in  References  K  and  Q  and  which  agrees  with  experimental  results. 

The  use  of  the  constant  property  relations  in  connection  with  the 
reference  enthalpy  method,  which  was  developed  for  the  laminar  boundary  layer, 
also  can  be  applied  for  the  turbulent  boundary  layer. 


Schultz-Grunow’s  equation  is  used  for  the  skin  friction  relationship 
for  turbulent  flow,  as  follows: 


0.370 
(tog10  ig 


2.584 

e 


(A-27) 


The  recovery  factor  is  represented  by  an  expression  similar  to  that 
for  laminar  flow 

r  =nv 


A- 7 


(A-28) 
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The  determination  of  the  heat  transfer  coefficient  still  is  based  upon  equation 
(A-lS)  and  upon  h 

<A‘29) 

where: 

hj^  =  h  Cp  (based  upon  enthalpy) 

All  property  values  must  be  introduced  into  equation  (A-27)  and  (A-28)at  the 
reference  enthalpy  (A-25),  and  p  in  equation  (A-29)  at  recovery  temperature. 

The  recovery  factor  now  is  based  upon  enthalpies, 

i  -  i 

=  -Vs-0-  2  g  J  (A-30) 


and  the  final  equation  for  the  temperature  change  of  the  skin  reads 

-$■  *  _ (o.i8a){jsi(.vjj _  «  .  1  +  *  vo  )  -  <  a  V 

4t  W(%vF My)^  ^  B5S5  — a- 

This  equation  can  be  solved  with  any  of  the  known  integration  methods. 


(A-31) 


g.  ‘  Stagnation  Point  Heat  Transfer  --  Generally,  missiles  designed  to 
fly  at  high  velocities  do  not  have  pointed  nose  cones  nor  sharp  leading  edges, 
which  would  have  extremely  high  local  heat  transfer  rates.  Instead,  noses  and 
leading  edges  are  rounded.  Therefore,  the  two  most  common  forms  appearing  in 
missile  design  (cylindrical  and  spherical  surfaces)  are  discussed. 

(l)  Laminar  Flow  in  Stagnation  Region  —  The  expressions  describ-  ( 
ing  heat  transfer  at  the  stagnation  point  in  laminar  flow  are 
described  in  Reference  R.  For  the  cylindrical  surface, 


=  n.^7  f-A-V3*8  x  D  1 

'  V  V  ^ r  7RT 


|  0  x  p>)  Stac 
'p*  ^ /a> 


(A-32) 


and  for  the  spherical  nose 


St.  -  0.736  (-£*•  ^5  ^  ta,  (A-33) 

The  property  values  for  equations  (A-32)  and  (A-33)  are  introduced  into  the 

(a)  Prandtl  number  at  stagnation  temperature. 

(b)  Reynolds  number  at  free  stream  temperature. 


The  equation  for  the  heat  transferred  to  the  wall  then  is 

q.  =  ST  x  x  V  (i  .  -  i  )  -  e  c  T  4 

eo  rct)  «  stag  w  w 

TO5 

The  term  V  V«  is  described  as 

HTTP  4  J  a  x  p~ 

^  Vc  ^  P  stag 

and  since  across  a  normal  shock 


(A  *3*0 


(A-35) 
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Poo 

it  follows  that 


p  stag  (K  +  1)  x  pi  +  K-l  (K-l)  +  2 

"  =  (K  -  l)lt  =  2  L  “  2  K  Tig  -  (k-iT 


1 


_D  =  %  6  C(K-l)  +  2] 

p  K-l  (K-l)  *4  +  2  ~| 

/  (K+l)  & 

L  1  +  2  2K^  -  (K-l)  J 

(A-3  6) 


(A-37) 


Experimental  verification  for  the  validity  of  this  equation  may  "be  found  in 
Reference  E. 


Since  the  Reynolds  numbers  of  the  local  flow  in  the  vicinity  of  the 
stagnation  point  are  low,  it  is  expected  that  the  flow  in  this  region  is 
laminar.  In  that  case,  it  follows  that  the  maximum  heat  transfer  occurs  at  the 
stagnation  point. 


(2)  Turbulent  Flow  in  the  Stagnation  Region  —  With  increasing 
distance  from  the  stagnation  point,  it  is  possible  that  the  flow  becomes 
turbulent.  For  this  condition,  the  following  relationship  is  taken  from 
Reference  R. 


(A-38) 


A  is  a  constant  having  the  following  values: 

A  *  O.CA-O  for  cylindrical  surfaces 

A  =  O.Olj-2  for  spherical  surfaces,  I)  =  diameter  of  nose 

The  rate  of  heat  transfer  to  the  wall  is: 

<1*  -  St.  •  poo  •  V«  (ir  -  i^)  -  e  •  a  •  T^a  (A-39) 

Property  values  are  introduced  as  follows: 

(a)  Reynolds  number  at  free  stream  conditions. 

(b)  Prandtl  number  at  conditions  existing  at  the  outer  edge  of  the 

foundary  layer.  ^ 


The  skin  temperature  can  be  calculated  in  the  same  way  as  in  preced¬ 
ing  paragraphs. 
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APPENDIX  B 

AERODYNAMIC  HEATING  TEST  FACILITIES 


The  three  "basic  inodes  by  which  heat  can  be  transferred  from  one 
medium  to  the  other  are  as  follows: 

Conduction:  The  heat  is  transferred  as  internal  energy  from  one 
molecule  to  another.  It  can  occur  in  solid  bodies,  fluids,  or  gases. 

Convection:  Motions  of  a  macroscopic  nature  may  exist  in  fluids  or 
gases  and  the  heat  is  carried  along  with  the  motion  in  form  of  internal  energy. 

Radiation:  Solid  bodies,  as  well  as  fluids  and  gases,  can  radiate  and 
absorb  thermal  energy  which  is  transmitted  in  form  of  electromagnetic  waves. 

In  technical  processes,  all  three  modes  of  heat  transfer  usually 
occur  simultaneously  with  varying  degrees  of  participation  for  each  particular 
mode. 


The  methods  of  heat  transfer  used  in  high  temperature  test  facilities 
can  be  grouped  into  two  broad  categories;  nonconvective  and  convective  heat 
transfer.  Reference  V  presents  a  detailed  description  of  the  various  methods 
and  facilities  available  in  these  categories,  which  will  be  discussed  briefly 
in  the  following  sections. 

/  \  a.  Nonconvective  Heating  Facilities  —  The  application  of  noncon- 

vective  heating  equipment  ranges  from  small  items  being'  tested  for  material 
properties  to  tests  of  large  structural  specimens  and  complete  flight  vehicles. 
The  purpose  of  the  investigation  is  to  study  the  effect  of  high  temperature  on 
the  degradation  of  material  properties,  creep,  fatigue,  and  thermal  stresses. 
These  investigations  can  be  performed  under  a  variety  of  test  conditions,  such 
as  heating  without  application  of  external  loads,  combined  heating  and  static 
loading  for  short  and  long  periods,  combined  cyclic  heating  and  loading,  and 
programmed  heating  and  loading  according  to  the  mission  profile.  Table  B-l 
illustrates  the  capabilities  of  some  of  the  existing  nonconvection  heating 
facilities. 

The  most  widely  used  method  of  heat  transfer  today  is  the  radiant 
heating  method.  It  resembles  the  actual  process  of  aerodynamic  hepJ  ig  quite 
well,  since  it  also  is  a  surface  heating  process.  A  typical  arrangement  is 
shown  in  Figure  B-l.  From  the  consideration  of  the  mathematical  description 
of  the  physical  process  two  approaches  for  rapid  changes  of  heat  flux  have  been 
established.  These  are:  (a)  To  heat  certain  elements  with  a  high  thermal 
inertia  to  a  constant  temperature,  (b)  to  obtain  changes  in  the  heat  flux  by 
programming  the  distance  between  the  test  item  and  the  heating  elements.  A 
typical  heater  of  this  category  is  the  one  developed  at  MIT.  It  employs  sill- 
con-carbon  rods  ("Globar"-rods)  which  are  mounted  on  a  movable  wall  of  a  closed 
chamber  and  are  re si stance -heated  to  approximately  2800°F  with  the  chamber 
closed.  The  wall  of  the  chamber  is  then  swung  out  to  face  the  test  specimen 
and  the  desired  temperature  profile  on  the  test  specimen  is  obtained  by  program- 
ming  the  distance  between  test  specimen  and  radiator,  (c)  To  keep  the  distance 


B-l 
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between  the  test  specimen  and  the  radiator  constant,  (d)  To  obtain  rapid 
changes  in  heat  flux  by  using  heating  elements  with  low  thermal  inertia  and 
programming  the  power  input  to  the  elements.  A  typical  element  of  this  kind 
is  the  quartz  tube  heat  lamp.  Because  of  its  versatility  and  thermal  qualities, 
it  is  used  extensively  in  the  laboratory  for  simulation  of  the  thermal  heat 
flux  of  aerodynamic  heating. 

Table  B-l  Capabilities  of  Nonconvective  Heating  Facilities 


Equipment 

Specimen 

Size  (*) 

Heating 

Rate 

r  bto  1 

Temperature  . 

Rate 

p  deg.  F  ■, 

Maximum 

Temperature 

[deg.  F] 

Lsec/fta  J 

^  sec  .  J 

Radiant  heaters 
Quartz -lamp 
Radiator 

large 

150 

450  (**) 

3,200 

Arc  Image 

2,600 

% 

> — ' 

0 

0 

0 

CO 

Furnace 

small 

7,000 

Furnaces 

large 

- 

- 

2,000 

small 

- 

5,000 

Electrical-Re- 

si stance  heaters 

small 

350 

1,500 

(***) 

Induction 

large 

3,000 

9,000  (**) 

(***} 

Heaters 

small 

25,000 

75,000  (***) 

(**#) 

Electron  beam 

heaters 

small 

>100,000 

>300,000 

(#**) 

Notes: 


(*)  Large  specimens  —  typically  greater  than  one  square  foot  of 
heated  surface. 

(**•)  Uniform  heating  of  0.1- inch-thick  aluminum. 

(***)  Maximum  temperature  in  excess  of  melting  temperature  of  material. 

One  of  the  most  commonly  used  lamps  consists  of  a  coiled  tungsten 
wire  filament  in  an  argon  atmosphere,  which  is  enclosed  by  a  quartz  envelope  of 
0.375  inch  diameter.  The  lamps  can  be  assembled  in  arrays  of  different  config¬ 
uration  to  fit  the  shape  of  a  particular  test  specimen.  For  maximum  thermal 
flux,  the  lamps  are  closely  spaced  in  each  of  two  rows,  where  one  row  is  placed 
behind  the  other  in  such  a  manner  that  the  filaments  of  the  back  row  lamps  fill 
the  spaces  between  the  filaments  of  the  front  row  lamps.  Heating  rates  that  a 

are  attainable  with  quartz  lamp  heat  radiators  are  about  150  [ — — T]  with 

sec.  ft  ( 
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cylindrical  radiators. 

The  lamps  are  available  in  various  lengths  from  5  to  50  inches  with 
two  types  of  filaments  with  power  ratings  of  100  and  200  vatts  per  inch  at  the 
rat  id  voltage.  They  can  also  be  ope.-ated  on  a  short  time  basis  at  100  percent 
overvoltage  and  then  deliver  30C  and  bOO  watt3  per  inch,  respectively. 

The  thermal  inertia  of  the  quartz  tube  lamps  is  small,  therefore,  the 
thermal  flux  can  be  varied  rapidly  and  accurately.  The  type  of  radiator 
control  widely  in  usd  today  is  the  ignitron,  which  has  fast  response  and  which 
is  adaptable  to  a  wide  range  of  electrical  loads. 

The  thermal  capabilities  of  cooled  and  uncooled  quartz  lamp  heat 
radiators  are  compared  with  the  thermal  simulation  requirements  of  ICBM  and 
manned  reentry  vehicles  in  Figure  B-2,  which  also  were  taken  from  Reference  V. 
As  shown,  quartz  lamps  can  generate  the  thermal  requirements  for  a  large  por¬ 
tion  of  the  manned  vehicle  environment,  but  are  inadequate  for  ICBM  simulation. 
For  this  portion  of  the  diagram,  which  required  heat  transfer  rates  substan¬ 
tially  greater  than  those  available  from  quartz  lamps,  the  arc -image  or  solar 
furnace  is  used. 

In  the  arc-image  furnace,  radiation  from  the  hot  gases  in  the  anode 
crater  of  a  carbon  arc  is  utilized  to  achieve  heating  rates  of  several  thousand 
BTU/fts sec .  The  application  of  arc -image  furnaces  is  limited,  however,  because 
the  area  which  can  be  irradiated  is  too  small  for  the  purpose  of  thermal  test¬ 
ing  of  aerodynamic  configurations.  The  radiant  heating  rates  can  be  increased 
by  operating  the  carbon  arc  in  an  elevated  pressure  atmosphere,  and  fluxes  of 
2600  BTU/ftr  sec  have  been  reported.  Reference  W  provides  more  detailed  inform¬ 
ation. 


The  greatest  usefulness  fur  conventional  furnaces,  including  arc- 
imaje  furnaces  is  in  the  field  of  materials  research.  They  are  not  suitable 
for  applications  in  which  temperature  gradients  are  important.  Lower  tempera¬ 
ture  furnaces  (to  2000  degrees  Fahrenheit)  usually  can  accommodate  small  test 
items,  while  higher  temperature  furcaces  generally  are  of  a  tubular  design  and 
restrict  specimen  size  and  shape. 

In  the  category  of  electrical  resistance  heaters,  it  is  ’..dcessary  to 
distinguish  between  two  types  of  heating  equipment;  one  which  directly  heats 
the  specimen  by  utilizing  the  specimen's  electrical  resistance,  and  one  which 
indirectly  heats  the  test  item  through  the  use  of  a  resistance  heated  blanket 
in  direct  contact  with  the  specimen.  The  latter  method  is  obsolete,  whereas 
the  direct  resistance  heating  technique  is  widely  used  for  testing  metallic 
materials. 


The  main  advantage  of  induction  heating  is  its  ability  to  develop 
high  heating  rates  and  temperatures  in  the  test  item.  However,  in  contrast  to 
radiant  heating,  induction  heating  is  not  a  surface  phenomenon.  Depending  upon 
the  frequency  of  the  heating  current  and  the  resistivity  of  the  material,  the 
heating  is  concentrated  at  a  certain  depth  below  the  surface.  In  addition,  the 
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FIGURE  B-2.  THERMAL  SIMULATION  CAPABILITY  OF  TWO  TYPES  OF  RADIANT  HEATERS 
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practical  application  of  induction  heating  for  aerodynamic  heating  simulation 
is  limited  by  the  physical  arrangement  of  the  heating  coils. 

Heat  transfer  rates  in  excess  of  100,000  BTU/fts  sec  have  been  ob¬ 
tained  over  small  ureas  by  bombarding  the  surface  with  a  direct  beam  of 
electrons.  The  technique  is  relatively  new  and  its  full  usefulness  in  the 
simulation  of  aerodynamic  heating  has  not  been  determined. 

b.  Convective  Heating  Facilities  —  Convective  heating  bears  the 
closest  resemblance  to  actual  aerodynamic  heating.  It  is  possible  to  achieve 
duplication  of  true  velocity,  enthalpy,  and  heat  transfer  rates  in  large 
facilities  for  simulation  of  low  supersonic  velocity  flight.  In  the  region  of 
high  supersonic  and  hypersonic  flight  simulation,  the  size  of  the  test  specimen 
and  the  duplication  of  parameters  are  restricted.  Table  II  illustrates  some  of 
the  existing  convective  heating  facilities. 


Table  II.  Capabilities  of  Convective  Heating  Facilities 


Facility 

Specimen 
Size  (*) 

Heat  Transfer 
Rate 

r  BTU  , 

Stream 

Enthalpy 

r  BTU 

Total 

Temperature- 

[deg.  F] 

^sec/ft"  ■*(##) 

[  lb  J 

Resistance  heated 
Tunnels 

small 

55 

550 

2,200 

Ceramic  heated 
Tunnels 

small 

300 

1,100 

4,000 

Combustion  heated 
Tunnels 

large 

180 

1,400 

4,000 

Electric  heated 
arc  facilities 

small 

1,000 

18,000 

16,000 

Notes: 

(*)  large  Specimen-typically  greater  than  one  foot  in  diameter. 
(**)  To  models  of  one-inch  hemispherical  radius. 


For  a  moderate  range  of  velocities,  the  energy  in  the  air  stream  of  a 
supersonic  tunnel  can  be  supplied  from  resistance  heated  heat  exchangers.  The 
maximum  operating  temperature  of  materials  of  the  heat  exchangers  limits  the 
range  of  operation  to  2000  degrees  Fahrenheit  for  stainless  steel  and  2500 
degrees  Fahrenheit  for  Kanthal  (Swedish  iron  alloy).  A  heat  exchanger  of  the 
latter  type  exists  at  Langley  Research  Center  in  a  continuous  flow  Mach  12 
tunnel  and  provides  a  stream  with  a  total  temperature  of  2,250  degrees 
Fahrenheit.  (  ) 
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In  a  current  development  effort,  porous  carbon  is  heated  to  5000 
degrees  Fahrenheit  and  nitrogen  is  passed  through  it  at  low  velocities  to 
permit  heat  absorption  without  erosion  of  the  carbon.  The  nitorgen  may  be 
used  as  a  test  medium  or  oxygen,  heated  separately  in  a  metallic  heat  exchanger, 
may  be  added  downstream  of  the  carbon  heat  exchanger  to  provide  synthetic  air 
at  temperatures  substantially  higher  than  can  be  obtained  from  metallic  heat 
exchangers  alone. 


Ceramic  heated  tunnels  provide  another  technique  for  obtaining  high 
temperatures  in  wind  tunnel  air  streams.  The  ceramic  bed  consists  of  pebbles 
of  aluminum  oxide,  zirconia,  or  other  ceramic  high  temperature  materials. 

High  temperature  combustion  gases  are  passed  through  the  bed  until  the  pebbles 
attain  the  desired  temperature.  Compressed  air  then  is  blown  through  the 
ceramic  bed,  heated  to  high  temperature,  and  expanded  through  a  nozzle  to  high 
supersonic  velocities.  Temperatures  achieved  with  a  large  version  of  the  cera¬ 
mic  heated  air  jet  were  approximately  4000  degrees  Fahrenheit.  A  serious  draw¬ 
back  of  ceramic  heated  facilities  becomes  evident  after  some  period  of 
operation.  The  cyclic  heating  fractures  the  pebbles  and  the  high  temperature 
air  stream  becomes  contaminated  with  ceramic  dust  particles,  which  in  turn, 
erodes  the  test  models. 


Combustion  provides  an  economical  means  of  heat  addition  to  a  gas 
stream.  A  slight  deficiency,  however,  is  the  presence  of  high  temperature 
combustion  products  instead  of  high  temperature  air.  This  can  be  circumvented 
by  using  fuels  which  yield  products  of  combustion  with  properties  approximating 
those  of  air.  An  advanced  type  of  facility  employing  combustion  heating  is  the 
high  temperature  structures  tunnel  at  the  Iangley  Research  Center.  Natural 
gas  or  ethylene  is  burned  in  a  high  pressure  combustor  and  the  combustion 
products  are  expanded  to  a  velocity  of  Mach  7  with  stagnation  temperatures  to 
4,000  degrees  Fahrenheit.  The  stream  of  combustion  products  contains  79  per¬ 
cent  Nitrogen,  water  vapor  up  to  10  percent,  oxygen,  and  carbon  dioxide.  The 
power  in  the  stream  is  approximately  one  gigawatt  (lO^  watts). 


Another  facility  at  the  Langley  Research  Center  is  a  9-  by  l6-foot 
tunnel  which  employs  combustion  of  a  fuel  as  the  source  of  heat  energy,  but  the 
energy  is  transported  to  and  stored  in  a  stainless  steel  heat  sink.  Compressed 
air  is  then  passed  through  the  sink  and  heated  to  600  degrees  Fahrenheit. 
Although  this  is  a  moderate  temperature,  the  tunnel  is  capable  of  developing 
true  aerodynamic  heating  and  loading  of  Mach  3  speeds  at  35>000  feet  altitude 
on  relatively  large  models. 


A  major  breakthrough  in  the  effort  to  increase  the  energy  in  air 
streams  of  ground  facilities  was  the  use  of  electric  arcs  to  heat  the  test 
medium.  Although  a  great  extension  of  the  simulation  capabilities  of  wind 
tunnels  was  achieved,  arc  heated  facilities  as  yet  do  not  provide  complete 
simulation  over  the  entire  flight  region.  Because  of  disassociation  and 
ionization  in  high  temperature  air,  the  pressures  required  in  the  arc  chamber 
to  produce  gas  streams  of  reentry,  Mach  numbers  become  impracticably  high.  As 
a  result,  arc  heated  chambers  can  produce  convective  heat  transfer  rates  of 
over  1,000  BTU/fta  sec  and  enthalpies  of  18,000  BTU/lb  in  air  and  can  thus 
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simulate  these  parameters  over  most  of  the  flight  regions  for  manned  and  un¬ 
manned  vehicles ,  although  flight  Mach  number  simulation  must  be  omitted. 

To  provide  simulation  at  positions  on  reentry  bodies  other  than  the 
stagnation  region,  more  sophisticated  techniques  have  been  developed  in  which 
the  manipulation  of  subsonic  or  low  supersonic  streams  provides  laminar  or 
turbulent  flows  and  simulates  changes  in  heat  transfer  rate  and  flow  velocity 
over  extended  regions  away  from  the  stagnation  point. 

c.  Recent  Developments  —  Two  recent  developments  in  the  areas  of 
radiative  and  convective  heating  offer  promise  of  high  gas  velocities  and  high 
radiant  energy  fluxes.  They  are  the  magneto-gasdynamlc  accelerator  and  the 
laser,  which  are  still  in  their  initial  development  stages  and  not  yet  suitable 
for  high-power  application. 

The  operating  principle  of  the  magneto-gasdynamic  heater  is  to  accel¬ 
erate  ionized  gases,  such  as  those  produced  by  a  conventional  arc-heater,  by  a 
body  force  applied  to  the  charged  particles  in  the  stream.  As  a  result,  the 
gas  is  accelerated  in  a  given  direction  without  a  significant  increase  in 
static  stream  temperature.  Theoretically,  this  device  could  accelerate  the  gas 
to  any  desired  velocity;  however,  experiments  show  that  it  is  difficult  to 
increase  the  stream  velocity  by  orders  of  magnitude  without  employing  tremen¬ 
dous  amounts  of  energy. 

The  laser  is  a  source  of  nonthermal  radiation  in  which  most  of  the 
energy  is  emitted  as  a  coherent  beam  of  radiant  energy  which  can  be  focused  on 
a  spot  of  extremely  small  size.  Radiant  energy  flux  densities  in  the  spot 
surpass  everything  that  has  been  obtained  so  far  in  the  laboratory. 
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AFEENDDC  C 
INSTRUMENTATION 


A.  High  Temperature  Measurements  —  Temperature  measurements  on  a 
metallic  test  specimen  can  "be  performed  with  a  reasonable  degree  of  accuracy 
by  thermocouples  up  to  approximately  2000  degrees  Fahrenheit.  Above  this 
limit,  thermocouple  measurements  become  more  difficult  and  less  accurate. 
Reference  U  presents  a  detailed  study  of  errors  in  surface  temperature 
measurements  made  with  thermocouples,  when  both  the  thermocouple  wires  and  the 
surface  are  exposed  to  radiant  energy.  The  dependency  of  accurate  results  upon 
various  parameters  is  discussed  and  recommendations  are  made  for  the  proper  use 
and  installation  of  thermocouples. 

For  higher  temperatures,  various  types  of  pyrometers  are  suitable  for 
measurement  of  a  steady-state  environment,  but  become  less  accurate  as  the 
temperature  rate  increases. 

B.  Strain  Measurements  at  High  Temperatures  —  One  of  the  important 
aspects  of  the  structural  behavior  of  a  test  specimen  is  the  measurement  of 
strain.  Techniques  for  conventional  strain  gages  (installation  and  measure¬ 
ment)  are  well  established  for  a  normal  temperature  environment,  but  become 
exceedingly  difficult  at  temperatures  above  a  few  hundred  degrees  Fahrenheit. 

At  these  temperatures,  special  techniques  are  employed  to  measure  strain;  for 
instance,  optical  measurements  of  displacement  of  reference  marks  on  the  heated 
surface,  off-surface  strain  gages,  measurements  on  the  cold  side  (if  available) 
of  a  test  specimen,  or  the  half -bridge  strain  gage  described  in  detail  in 
Reference  X.  The  principle  underlying  the  operation  of  the  half-bridge  strain 
gage  is  to  induce  theoretically  equal  temperature  effects  in  each  component  of 
a  half-bridge  which  cancel  when  the  elements  are  connected  in  a  standard  wheat- 
stone  bridge  circuit.  Solutions  are  developed  for  the  compensation  for  effects 
such  as  thermal  coefficient  of  resistance  of  the  strain-sensitive  wire,  the 
self-heating  of  the  gage  caused  by  the  electrical  current  flowing,  the  strain 
sensitivity  of  each  arm,  leakage  resistance  of  the  cement,  and  the  lead 
resistance . 

C.  Gas  Property  Measurements  in  Convective  Heating  Facilities  --  The 
use  of  convective  heating  facilities  requires  the  measurement  of  additional 
parameters  related  to  the  properties  of  the  gas  stream.  These  parameters  are 
the  enthalpy,  velocity,  and  pressure  of  the  gas  and  the  heat  transfer  rate 
experienced  by  the  test  specimen.  They  are  discussed  under  "Measurements  in 
Arc-heated  Air  Streams"  in  Reference  V.  According  to  this  reference,  pressures 
easily  can  be  measured  using  standard  pressure  transducers.  Probes  used  for 
dynamic  pressure  measurements  must  be  water  cooled.  Velocities  can  be  measured 
directly  or  can  be  calculated  from  pressure  and  temperature  data.  Mach  numbers 
can  be  obtained  by  reading  photographs  of  the  shock  waves  formed  on  wedges 
immersed  in  the  gas  stream. 

Several  methods  have  been  used  successfully  for  obtaining  the  surface 
heat  transfer  rate  which  is  one  of  the  most  important  parameters  in  high  temper- 
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ature  testing.  The  simplest  form  of  such  a  measuring  device  is  a  thin  metallic 
disk  with  a  thermocouple  attached  to  the  center.  When  the  disk  is  immersed  in 
the  gas  stream,  the  rate  of  temperature  rise  is  a  direct  measure  of  the  heat 
transfer  rate. 

For  the  measurement  of  stream  enthalpies  and  temperatures,  several 
direct  and  indirect  methods  have  "been  developed.  The  gross  enthalpy  may  be 
determined  by  measuring  the  total  flow  weight  and  total  energy  of  the  stream. 

The  latter  can  be  calculated  by  subtracting  the  energy  lost  in  the  arc  appara¬ 
tus  from  the  total  power  input  to  the  arc  heater.  Another  way  to  determine 
the  total  energy  is  to  collect  the  entire  gas  stream  in  a  tube  and  measure  the 
energy  absorbed  by  the  tube  in  cooling  the  gas.  A  useful,  indirect  method 
for  determining  the  enthalpy  or  temperature  as  a  function  of  position  in  the 
stream  is  based  on  the  heat  transfer  rate  theory  of  Fay  and  Riddell.  Stagna¬ 
tion  point  heat  transfer  rate  on  a  body  is  calculated  as  a  function  of  stagna¬ 
tion  enthalpy  and  density,  and  local  velocity  gradient.  Newtonian  flow  theory 
is  used  to  calculate  the  velocity  gradient  in  supersonic  flow,  whereas  in  the 
subsonic  region,  potential  flow  theory  is  applied.  Results  are  meaningful 
only  for  Reynolds  numbers  above  1000. 

Another  direct  measurement  of  enthalpy  at  particular  locations  in  the 
stream  can  be  made  by  using  a  small  gas  sampling  probe  which  draws  off  a 
measured  weight  flow.  The  enthalpy  is  determined  from  the  heat  extracted  in 
the  probe. 

Spectrophotometry  is  one  of  the  methods  used  in  the  direct  measure¬ 
ment  of  stream  temperatures.  It  utilizes  the  spectral  intensity  distribution 
of  energy  radiated  by  the  stream  and  it  generally  yields  temperatures  indica¬ 
tive  of  one  degree  of  freedom  of  the  gas.  These  temperatures  describe  the 
energy  state  of  a.  .1  components  in  the  gas,  provided  the  gas  is  in  thermodynamic 
equilibrium. 

D.  Data  Handling  Equipment  —  Si  dynamic  test  applications,  the  rapid 

collection,  reduction,  and  interpretation  of  data  is  essential  to  the  successful 
completion  of  the  program.  This  can  be  accomplished  only  through  the  use  of 
automatic  electronic  data  handling  equipment.  Various  types  of  automatic  data 
processing  equipment  are  available,  but  it  would  exceed  the  scope  of  this 
section  to  discuss  all  of  them  in  detail. 

The  important  aspects  for  the  acquisition  of  a  data  handling  system 
as  presented  in  Reference  T  are  listed  as  follows: 

a.  Sufficient  ciiannel  capacity  and  sampling  speed,  in  the  commutated 
mode,  to  support  test  requirements  ranging  from  static  loading  with  steady- 
state  temperatures  to  dynamic  transient  loadings  and  temperature  with  the 
option  of  sacrificing  channel  capacity  for  increased  sampling  rate. 

b.  The  capability  of  programming  the  system  to  operate  in  a  fully 
automatic  mode  when  events  are  occurring  so  rapidly  that  human  perception  or 
judgement  cannot  be  applied. 

c.  System  accuracy  such  that  the  transducers  contribute  the  major 
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sources  of  error. 

d.  Data  editing  capabilities  so  that  data  can  be  examined  for  its 
informational  content  prior  to  machine  processing. 

e.  Modularization  so  that  a  part  of  the  system  may  be  used  to 
obtain  data  from  one  test  article  concurrently  with  test  setup  operation  of 
the  remainder  of  the  data  system  on  a  second  or  third  article. 

f.  Isolation  capabilities  such  that  transducer  failure  or  malfunc¬ 
tion  on  one  channel  will  not  affect  the  operation  or  accuracy  of  any  other 
channel . 

g.  The  capability  of  processing  and  putting  out  a  large  volume  of 
test  data  within  a  short  period  of  time  for  examination  and  analysis  by  instru 
mentation  engineers. 

h.  Sufficient  growth  capacity  such  that  the  system  will  not  become 
obsolete  within  a  reasonable  time  period. 

i.  Component  or  module  accessibility  for  setup  and  maintenance 

purposes. 


To  show  the  capability  of  typical  system,  the  description  of  the  WADD 
elevated  temperature  and  fatigue  test  facility  data  handlihg  system  is  quoted 
from  Reference  T:  "This  system  collects  (by  commutation  and  with  a  variable 
sampling  rate),  records,  analyzes,  reduces,  displays,  and  stores  the  outputs  of 
various  sensing  devices  such  as  strain  gages,  load  cells,  deflection  potent io- ' 
meters,  and  thermocouples.  Electrical  signals  from  the  transducers  are  detected, 
conditioned,  and  transmitted  by  shielded  cables  to  a  data  control  and  processing 
area.  Here  the  data  are  visually  monitored,  recorded  on  magnetic  tape,  process¬ 
ed  and  stored  by  means  of  a  large  scale  digital  computer  and  a  high-speed 
printer-plotter.  Also  included  are  'limit -alarm'  and  'quick-look'  display 
features. 


"The  system  consists  of  eight  transmitter-multiplexer  units,  two 
monitor  consoles,  four  one-inch  tape  recorders,  four  one -half  inch  tape  record¬ 
ers,  one  system  control  console,  a  large  scale  digital  computer  and  a  high¬ 
speed  electronic  printer  plotter. 

"Total  channel  capacity,  for  the  full  scale  system,  is  1920  channels. 
A  basic  block  diagram  of  the  full  scale  system  is  shown  in  Figure  C-l.  The 
channel  capacity  can  be  increased  by  simply  adding  transmitter-multiplexer  and 
tape  recorder  units." 
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Figure  C-l  Data  Handling  System  Block  Diagram 


MTP  5-2-509 
2b  July  1967 


APPENDIX  D 

BASIC  TECHNIQUES  FOR  PROGRAMMING  MISSION  PROFILES 

The  basis  for  programming  thermal  profiles  is  equation  (A-8),  which 
is  solved  by  a  heat  control  computer.  The  principle  of  operation  is  demon¬ 
strated  through  the  use  of  a  radiant  heating  test  facility,  although  basically 
it  applies  to  different  types  of  heating  equipment. 

a.  Heat  Control  Computer  Purpose  —  The  purpose  of  the  heat  control 
computer  in  conjunction  with  electrical  power  control  equipment  (ignitrons)  is 
to  regulate  the  heat  applied  to  test  structures  to  duplicate  the  effect  of 
aerodynamic  heating  resulting  from  flight  at  supersonic  speeds.  The  method  of 
operation  presented  here  is  taken  from  Reference  T,  which  describes  part  of  the 
radiant  heating  facility  at  Wright  Air  Development  Division  (WADD). 

The  calculation  procedures  in  Appendix  A  show  that  the  analytical 
determination  of  the  instantaneous  specimen  temperature  Ty  is  tedious.  As 
used  in  the  computer  facility,  it  is  an  actual  test  derivative  from  the  control 
thermocouple  for  the  test  area  under  consideration  rather  than  by  calculation. 
For  test  purposes,  the  variation  of  the  heat  transfer  coefficient  h  and  the 
recovery  temperature  Tr  with  respect  to  time  are  required.  The  remaining  items 
are  known  for  a  specific  test  article  and  the  given  test  conditions.  Equation 
(A-8)  can  be  written  in  a  different  form,  as  follows: 

Hocal  -  if*  -  h  (^  -  <)  -  B  V 

where:  B  =  e*cr  (a  constant  for  a  specific  test  article) 

Further,  the  thermal  input  can  be  related  to  the  power  input  of  the  laboratory 
heat  source  by 

Q  =  YxExIx  10"3 

where: 

Q  =  heat  produced  in  the  electrical  process 

Y  =  a  constant  relating  K W  to  BTU/sec 

E  =  load  voltage 

I  =  load  current 

Reference  T  states,  ’’The  heat  control  computer  solves  equation  (A-8) 
instantaneously  using  as  inputs  the  precalculated  values  for  the  recovery 
temperature  Tr  and  the  heat  transfer  coefficient  h  together  with  the  measured 
surface  temperature  of  the  test  structure  Tw.  The  generated  output  signal 
controls  the  heat  source  to  produce  the  proper  heat  flow  rate  Q." 

b.  Computer  Operation  —  A  functional  block  diagram  of  a  forty- 
channel  heat  control  computer  is  shown  in  Figure  D-l  (Reference  T).  Curves 
are  plotted  from  the  precalculated  input  data.  Three  hundred  points  for  each 
of  ten  recovery  temperature  Tr  curves  and  forty  heat  transfer  coefficient  h 
curves  are  punched  in  a  seven  place  code  on  paper  tape  and  transferred  in 
digital  form  by  means  of  tape  reader  and  write  circuitry  to  a  magnetic  storage 
drum.  Through  use  of  the  read  circuitry,  each  curve  can  be  plotted  on  a  strip 
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chart  recorder  to  verify  visually  the  accuracy  of  the  information  that  is 
stored  on  the  drum.  The  read  circuitry  scans  points  of  equal  abscissa  from  the 
curves  during  each  drum  revolution  as  programmed  by  a  time  base  selector. 
Although  the  time  base  range  covers  thirty  seconds  to  one  hour  in  increments  of 
10  seconds,  the  scanning  may  be  held  at  any  of  the  abscissas  up  to  eight  hours. 
The  digital  outputs  of  the  read  amplifiers  are  converted  to  analog  voltages  by 
means  of  D-A  converters.  The  forty  h  curves  are  applied,  through  the  convert¬ 
ers,  directly  to  a  high  speed  rotary  multiplexing  switch.  The  ten  Tr  curves 
are  applied,  through  the  converters,  to  forty  ten-position  switches.  The 
selector  arm  on  these  switches  is  in  turn  applied  to  the  multiplexing  switch, 
thereby  allowing  one  of  ten  possible  Tr  temperature  selections  to  be  used  for 
each  of  forty  test  areas. 

Temperature  input  is  derived  from  a  thermocouple  through  a  low  level 
thermocouple  amplifier  circuit.  In  addition  to  amplification,  the  thermo¬ 
couple  is  linearized  so  that  its  output  signal  is  proportional  to  temperature. 

A  thermocouple  test  circuit  is  available  to  test  the  operation  of  each  ther¬ 
mocouple  channel.  Thermocouple  failure  circuits,  connected  to  the  low  level 
amplifiers,  indicate  any  open  circuit  in  the  thermocouple  input  circuits. 

The  radiation  factor  B  =  «*cr  is  set  manually  for  each  channel. 

The  signals  which  represent  Tr  and  -  Tw  are  summed  in  amplifier 
number  1.  The  output  signal  of  this  anrplifier  is  modified  by  the  signal 
representing  h  in  the  multiplier.  The  multiplier  output  is  a  voltage  propor¬ 
tional  to  -h  (Tr  -  Tw)-  The  representative  signal  -  B*  *  Tw  is  applied  to  the 
input  of  the  function  generator  to  produce  a  signal  proportional  to  BT^.  The 
signals  -  h  (Tr  -  Tw)  and  BTy4  are  added  in  summing  amplifier  number  2  and 
the  output  represents  h  •  (Tr  -  T^)  -  BTy4  which  is  proportional  to  the  desired 
power  Q.  This  output  is  commutated  by  the  multiplexing  switch  to  the  output 
holding  circuits.  A  combined  area  and  heat  transfer  efficiency  factor,  K,  is 
introduced  and  this  output  is  compared  with  either  the  rate  of  temperature 
rise  ( c ♦  s * 6 •  *][~w  )  or  the  power  consumed  (Y*E* I)  feedback  loops  (depending  upon 
dt  . 

the  mode  of  operation)  in  summing  amplifier  number  3-  The  resulting  signal  is 
used  to  control  the  power  to  the  heat  source. 

In  the  c’s*£W  .  6,  or  derivative  feedback  mode,  a  signal  from  the 

'  dt" 

thermocouple  linearizer  is  differentiated.  This  constant  c*s*6  is  set  in  and 
the  signal  is  entered  in  amplifier  number  3*  In  the  power  feedback  mode, 
signal  voltages  representing  the  voltage  and  current  in  the  lamps  are  multi¬ 
plied  to  produce  a  signal  which  represents  the  power  used  by  the  lamps.  This 
signal  is  averaged  by  a  filter  network  to  represent  the  actual  power  per  cycle. 
The  output  of  this  network  is  modified  by  the  constant  conversion  Y  and  is  sent 
to  summing  amplifier  number  3- 

If  the  wall  temperature -time  history  of  the  test  items  has  been 
established,  then  the  h  curves  may  be  substituted  by  the  T„  curves  in  this  mode 
of  operation.  (See  Figure  D-2)  The  signals  representing  T^  are  fed  directly 
from  the  D-A  converter  to  the  output  holding  circuit.  T^  is  multiplied  by  K 
and  the  signal  KT  is  then  compared  with  the  actual  value  of  Tw  amplifier  and 
fed  through  c*s.fi.  The  generated  error  signal  controls  the  power  units. 
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FIGURE  0-1.  FUNCTIONAL  DIAGRAM  OF  WADD  (G.E.)  40  CHANNEL  HEAT  CONTROL  COMPUTER 


